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An  Instability  Theory  of  Ice-Air  Interaction  for  the  Formation 

of  Ice  Edge  Bands 

p.  C.  Chu 

Deportment  of  Oceanography,  Naval  Postgraduate  School,  Monterey.  California 

Icc  edge  bands  arc  well-observed  phenomena  at  the  ice  edge  of  the  world  ocean  at  times  when  winds 
are  blowing  off  the  ice.  Differing  from  the  wave  radiation  theorv-,  water  wave  theory,  wind  wave  theory, 
and  ice-water  coupling  theory,  we  propose  an  ice-air  feedback  mechanism  to  generate  the  ice  edge  bands. 

To  climmate  ocean  effects,  the  water  is  set  motionless.  Thermally  generated  surface  winds,  blowing  from 
ice  to  water  (ice  breezel  with  some  deflection  due  to  the  earth's  rotation,  force  the  drift  of  ice  floes  in  the 
marginal  ice  zone  iMIZl.  By  changing  the  surface  icmperalurc  gradient  the  ice  motion  feeds  back  on  the 
surface  winds  A  coupled  ice-air  model  similar  to  the  author's  earlier  paper  is  employed  to  discuss  the 
instability  propierties  of  such  a  feedback  mechanism.  The  linearized  governing  equations  are  then  solved 
as  an  eigenvalue  problem  with  two  model  parameters  mean  ice  thickness  H,  (0..^  m  <  H,  <  Z5  ml  and  a 
characteristic  surface  temperature  difference  over  icc  and  water  DT^  (4'C  <  DT^  <  20'  C).  The  model 
results  show  alternating  ice  divergence  and  convergence  zones  and  demonstrate  a  great  num'oer  of 
e-folding  acceleration  modes  exated  by  the  coupled  ice-air  system.  Ice  floes  are  transponed  from  the 
divergence  area  to  the  convergence  zone  and  are  eventually  rearranged  into  a  band  structure  The 
e-folding  acceleration  rate  of  ice  motion  depends  on  the  width  of  two  adjacent  convergence  zones  L  It 
increases  with  L'l  (where  L  is  the  length  of  twice  the  MIZ  width,  ix..  200  kml  for  10  km  <  /  <  200  km. 
and  then  remains  high  for  2  km  <  f  <  10  km  The  length  scale  (2-10  kml  of  the  high  acceleration  modes 
agrees  well  with  the  observed  ice  edge  bandwidth  (1-10  km).  The  e-folding  acceleration  rale  increases 
with  an  increase  in  DTf,  and  with  a  decrease  in 


1.  I-STRODUCnON 

Surface  bands  of  icc  floes  arc  generally  observed  at  icc  edges 
of  the  world  ocean  when  winds  arc  blowing  off  the  ice.  These 
bands  have  size  scales  of  I -10  km.  with  their  long  axes  orient¬ 
ed  approximately  normal  to  the  wind  direction.  The  icc  band 
spacing  varies  strongly  from  1  to  10  km.  Differing  from  wave 
radiation  theory  [H'adham.s.  1983],  water  wave  theory 
[.Muench  ei  aL  1983],  wind  wave  theory  [.Martin  et  al..  1983], 
and  icc-watcr  coupling  theory  [Hdkkinen,  1986].  wc  propose 
in  this  study  an  icc-atr  feedback  mechanism  for  the  generation 
of  ice  edge  bands  in  the  marginal  ice  zone  (MlZl  by  means  of 
an  icc-air  model  similar  to  that  treated  by  Chu  f  1986a]  except 
that  ( 1 )  no  horizontal  eddy  viscosity  v*  is  included  in  the 
boundary  layer  airflow  model  and  (2l  the  .solutions  arc  not 
restricted  to  the  first  Fourier  mode. 

MePhee  [1980]  mentioned  ice-ocean  and  ice-air  coupling  in 
a  review  paper.  A  possible  mechanism  for  the  meso.scale  ice-atr 
interaction  is  presented  in  Figure  I.  Surface  winds  generated 
by  the  surface  temperature  gradient  blow  from  ice  to  water 
(ICC  breeze!  with  .some  dcflecuon  due  to  the  earth's  rotation. 

However.  movement  of  ice  floes  in  the  .MIZ  changes  thermal 
conditions  near  the  surface  and  produces  a  diffcreni  surface 
temperature  gradient  across  the  ice  edge,  which  in  general  is 
less  than  or  greater  than  the  characteristic  temperature  gradi¬ 
ent  in  the  MIZ. 

The  ice-atr  interaction  model  depicted  in  the  subsec)uent 
sections  is  intended  to  determine  the  unstable  modes  of  ice 
motion  and  to  predict  the  length  scale  of  the  ice  band  spacing. 

2.  Thermally  Forced  Boi'ndary  Layer  Air  Floia- 

Overland  et  al.  [1983]  show  that  air  temperature  increases 
w'atcrward  in  the  Bering  Sea  MIZ  (Figure  2).  Such  a  surface 
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temperature  gradient  will  thermally  generate  local  atmospher¬ 
ic  flow  near  the  MIZ.  In  this  section  we  utilize  a  planetary 
boundary  layer  air  model  treated  by  Kuo  [1973]  and  Chu 
[1986a.  h.  c]  to  simulate  thermally  forced  boundary  layer  air¬ 
flow.  The  .X  axis  is  m  the  cro.ss-edge  direction,  and  the  y  axis  is 
parallel  to  the  ice  edge,  as  is  shown  in  Figure  3. 

The  potential  temperature  of  air  is  divided  into  two  parts:  a 
basic  state  sod  perturbation  0^‘.  The  basic  state  is  given 
by 

—  Oyff  -t-  (1) 

where  Ogo  is  the  basic  air  potential  temperature  at  the  surface 
and  jV  IS  the  Brunt-Vaisala  frequency.  The  subscript  asterisk 
means  dimensional  variables. 

As  discussed  by  Chu  [1986a].  the  waterward/iceward  migra¬ 
tion  of  the  MIZ  increases' decreases  the  surface  temperature 
gradient,  so  the  effects  of  ice  flow  on  the  surface  air  temper¬ 
ature  gradient  can  be  parameterized  as 

=  <nD7i,/L*)u,,(.x,.  ;,)  (2) 

where  DTo  is  the  characteristic  surface  temperature  difference 
across  the  MIZ  and  L  is  a  length  scale  twice  the  MIZ  width: 
subscripts  a  and  i  denote  air  and  ice.  respectively. 

The  coordinates  and  atmospheric  variables  are  nondimen- 
sionalized  by  setting 

(.v,.  z,.  r.l  =  {xL  :S.  tT]  s'  =  0,,'iQg^  =  (DTg/Ogois, 

<3) 

(u.,.  i-..,  u„.  p..!  =  Ulu^  v„  w.  6/L,  u,,  i;,) 

P.,  =  T^/0ao)P.  ^  =  ( »'./«) '  '• 

where  v  is  the  vertical  eddy  viscosity,  fi  is  the  angular  velocity 
of  the  earth's  rotation. 
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Fig.  I.  Ice-air-ocean  coupled  system. 


is  the  scale  of  icc  breeze  wind,  and  T  =  L/V  is  the  time  scale 
for  change  of  the  surface  temperature  gradient  due  to  move¬ 
ment  of  the  icc  edge.  If  we  assume  that  local  airflow  satisfies 
the  modified  Boussinesq  approximation  [ACuo.  1973]  then  the 
vorticity  equation,  the  momentum  equation  (both  in  the  y 
direction!,  and  the  heat  equation  for  air  disturbances  gener¬ 
ated  by  the  diflerential  surface  temperature  gradient  near  the 
MIZ  are  [Clm.  1986a] 

=fodvJcs  -  dsjdx  (5) 


T 


> 

Ice  Interior 

MIZ 

* 

* 

'  Open  Water 
* 

0 

Fig.  1  The  model  MIZ  and  the  coordinate  system. 

Ed^vjdz^  =  -f^^^|>J^z  (6) 

Ev^sJSz'  =  Ri  Btjtjdx  (7) 

In  the  foregoing. 

=  -Bil/Jdz  w,  =  B^Jdx 

(8) 

V*  =  S^/l}ci^/dx^  +  d^ldz^ 

where /o  =  sin  is  latitude,  and 

E  =  v/(2£li5^)  =  ^  Ri  =  (9) 

are  the  Ekman  and  Richardson  numbers,  respectively. 

The  model  variables  are  decomposed  into  Fourier  sine  or 
cosine  series: 


Fig  2.  Cross  sections  of  potential  temperature  perpendicular  to 
llic  ice  edge  in  the  Bering  Sea  for  the  three  ofl-icc  wind  ca.scs  [from 
Oivrland  er  at..  1983]. 


il/,{,x.  z,i)’=  Y,  I)  sin  knx  (10a) 

k>  I 

S.I.-0  =  Z  s-»(0  cos  (fO*) 

I*  t 

-ac 

“.Ir-O  =  Z  “«»(')  S'" 

1 

(lOc) 

'’ali  =  0  =  Z 

k*  I 

U,=  Y.  S*" 

(lOd) 

•Jt 

V,  =  ^  p^(f)  sin  kxx 

*•  I 

Introducing  a  new  variable  f,  representing  nondimensional  ice 
displacement  in  the  x  direction. 

LdQdi^u,  (II) 

expanding  it  into  Fourier  sine  series. 

c,(x- 0  =  Z  (12) 

t«  I 

and  substituting  (lOif)  and  (12)  into  (1 1),  we  obtain 

LdCJdt  =  Ui^  (13) 

Integrating  (2)  with  respect  to  x  and  t  after  nondimensionali- 
/ation.  we  get 

s»=-Kui(0  +  ^J/k  (14) 

where  fi^ik  =  1, 2,  •  •  • )  arc  the  integration  constants. 
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Eliminating  p,  and  from  (5H7).  we  find  that  the  stream 
function  satisfies  the  following  partial  differeniiai  equation: 

( l/4f?*,rs*  -f-  +  Ri  =  0  (15) 

We  solve  (15)  for  the  stream  function  and  obuiin  the  solu¬ 
tions  of  p,  and  Sj  from  (6)  and  (7)  after  substituting  4’,-  The 
local  airflow  is  thermally  forced  by  the  surface  temperature 
gradient,  as  is  indicated  in  (14).  Tbe  boundary  conditions  in 
the  vertical  direction  arc  derived  as  follows.  The  dependent 
variables  should  remain  finite  as  r  -♦  oo,  i-C-. 

lim  (|<&,|.  \c4>JB2\,  |pJ.  |.v^|)  <  x  (16) 

The  boundary  conditions  at  r  =  0  arc 
=  0  Siiijoz  =  M 

(17) 

<r> 

P,  =  M  cvjcz  s,  -  -  Z  +  W  COS 

*  =  I 

where  M  is  a  measure  of  the  effective  depth  of  the  constant 
stress  sublayer  [Kuo.  1973]:  Af  =  v/{C\,S),  where  C„  is  a  di¬ 
mensional  (meters  per  second)  air  drag  coefficient. 

Substituting  (lOw)  into  (15).  wc  obtain  the  following  sixth- 
order  ordinary  differential  equations  for  the  Fourier  coef¬ 
ficients 

4/o - Ri  <1/^=0  (18) 

The  general  .solutions  of  (18)  have  the  following  form: 

’)  =  (Si*  +  P»)  y,  exp  (a,^)  (19) 

where  the  eigenvalues  (j=  1.  2,  •  •  • .  6)  are  the  roots  of  the 
sixth-order  algebraic  equatioas: 

;.*■  -I-  4/oV.»  -  Ri  -  0  (20) 

According  to  the  upper  boundary  conditions  listed  in  (16),  we 
must  set  up  coefficients  that  corre.spond  to  those  eigenvalues 
with  the  positive  real  piarts  to  zero.  Consequently,  the  general 
solutioas  (19)  sati-sfying  the  top  boundary  conditioas  may  be 
written 


3 


0  =  rcu(')  +  w  E 


(21) 


where  the  eigenvalues  all  have  negative  real  parts.  Substi¬ 
tuting  (21 )  into  ( lOu).  we  obtain  the  stream  function 


z,r)=  y  a^jCCa  +  PJ  exp  sin  (22) 

Using  the  definition  of  a  stream  function  and  integrating  the 
momentum  equation  (6)  and  the  heat  equation  (7)  with  respect 
to  c  after  substituting  (22),  wc  find  that 


“.Uo  =  Z  "«*[sa  M  sin  fetex 
'’Jr-O  =  X  +  fiJ  S‘n 

k  =  i 

x> 

~  Z  /fj  20$  kxx 

l^l 


(23) 


where 


1  3 

=  Z  ^  -  ^0  Z 

j^i  /“I 


=  2knRi  Y 


(25) 


Substituting  solutions  (22)-(24)  into  the  surface  boundary  con¬ 
ditions  listed  in  (17),  we  get  the  following  four  algebraic  equa¬ 
tions  for  an  fc, : 


Z  =  0 


Y/.^/\  -  .VfA*j)a»,  =  0 
J»i 

■*  I  —  M} 

-%Z  +  =  0 

1 

3 

2knR,  Y,  =  -  I/* 


(26) 


1 


3.  Free  Ice  Drift  Model 

The  linearized  momentum  equations  for  a  free  drift  ice 
model  are 

3D 

=  -  C^Ui  +/Vi  -1  C  Z  “»»Ka(ri  +  Pti  sin  (22) 

■X> 

dvjdt  =  -C^^  -fu,  +  C.  Z  2»k[(^a(0  +  W  sin  (28) 

*- 1 

where  the  first  and  last  terms  on  the  right-hand  side  represent 
the  water  and  air  stresses  on  ice,  respectively.  Here 

=  P.CJUhHi  C  s  p^C^I{p,H,)  (29) 

C,.,  is  a  dimensional  (meters  per  second)  water  drag  coefficient 
(on  ice),  and  p„.  p,,  and  p„  arc  the  densities  of  air,  ice.  and 
water,  respectively. 

4.  Unstable  Modfs  of  Ice  Motion 
The  solutions  of  (27)-(28)  have  the  following  forms: 

V 

M/.X,  /)  =  X 

t 

3e 

r/.x,  t)  =  y  Pu  (®»ri  •'•in  kxx  (31) 

k*l 

Substituting  (30)  and  (31)  into  (27)  and  (28).  wc  obtain  the 
following  dispersion  relation: 

+  20i*  +  iCJ  -  Cji^  L 

-  (CjDlCji^  -f/bJ  =  0  (32) 

The  roots  e^{k  =  I,  2.  •  •  •)  of  the  cubic  equations  (32)  give  the 
c-folding  time  dependence  of  the  Itth  component  of  ice  veloci¬ 
ty.  The  standard  values  of  model  parameters  are  given  in 
Table  1.  The  instability  criterion  for  the  kth  mode  of  ice 
motion  in  the  MIZ  can  be  written 


(24) 


Re  (Oi)  <  0 


(33a) 
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TABLE  I.  Siandard  Model  Parameiers 


Parameter 

Value 

L 

200  km 

n 

0.7292  X  10-“*-' 

P. 

1.29  kgm“ 

c. 

3  X  10"^  m/s 

V 

5  mVs 

u 

9.81  ms’ 

p„ 

lO’  kg  m’ 

c.. 

1.76  X  I0~’  m.s 

<p 

65" 

270  K 

Pi 

910  kg  m’ 

the  fcth  mode  of  icc  velocity  decreases  with  time. 


neutral,  and 


Re  (ff^)  =  0 


(33b) 


Re  (o»)  >  0  (33c) 

the  lith  mode  of  icc  velocity  increases  with  time.  We  define  the 
time-increasing  mode  of  ice  velocity  as  an  unstable  mode. 

The  oscillation  criterion  for  the  /ith  mode  of  icc  motion  is 


given  by 

Nonoscillatory 

Im  (ffj  =  0 

(34fl) 

Oscillatory 

Im  (<7^)  9*  0 

(34b) 

We  compute  all  the  roots  of  (32)  for  different  values  of  the 
parameters  H,  and  OT^.  which  are  organized  into  two  groups; 

(1)  fixed  //j  (1.5  m).  DTo  =  4  C  8  C.  12X.  and  20'-C  and 

(2)  fi.xcd  DTo  (lO'C).  Ht  =  0.5.  I.  1.5.  2.  and  2.5  m.  We  obtain 

three  roots  <7,(k),  and  ffj(fc)  for  each  given  k.  and 

DTo-  Here  (7,(fe)  is  real  and  positive  for  all  given  fe,  Hj.  and 
DTo-  However,  c^ik)  and  ajik)  are  complex  with  negative  real 
parts.  Therefore  the  time-increasing  modes  of  icc  velocity  are 
those  eigenfunctions  associated  with  the  eigenvalues  a^^k): 


0  =  sin  knx 

c)  =  exp  r<7,(^)f]  sin  knx 


(35) 


which  indicate  that  the  ice  velocity  =  (iJ^j.  Sg)  has  an  e- 
folding  acceleration  and  a  sinusoidal-type  spatial  variation. 


Fig.  5.  The  a,^k)-k  diagram  for  =  1.5  m  and  (u)  DT^  =  4‘’C.  (b) 
DTo  =  rc.  (e)  DTo  =  I2'C  (d)  DTo  I6"C.  and  (el  DTo  =  ^O'C. 


5.  Model  Results 

Solutions  (35)  show  that  time  and  spatial  variations  of  ice 
velocity  arc  separable,  which  means  that  the  sinusoidal-type 
spatial  structure  does  not  change  with  time  and  that  the  e- 
folding  acceleration  docs  not  vary  with  space. 

Figure  4  indicates  the  spatial  variation  (sinusoidal)  of  the 
Ath  mode  icc  velocity  depicted  in  (35).  We  see  the  alter¬ 
nation  pattern  of  icc  convergence  and  divergence  zones.  The 
width  of  the  two  adjacent  convergence  zones  is  L/k.  Such  an 
ice  velocity  field  transports  ice  floes  from  the  divergence  area 
to  the  convergence  zones  and  eventually  forms  ice  edge  bands 
with  band  spacing  L/k. 

Solutions  (35)  also  show  that  this  convergenoe-divergcncc 
alternation  pattern  has  an  e-folding  time  variation.  The  e- 
folding  acceleration  factor  o-ifA)  denotes  how  fast  the  ice  veloc¬ 
ity  increases.  Figures  5  and  6  arc  the  o', -A  diagram  for  groups 
1  and  2  of  parameters  H,  and  DTf,,  respectively.  “T-type" 
curves  demonstrate  the  following  facts. 

1.  For  any  given  parameters  //,  and  D7J,,  the  e-folding 
acceleration  factor  o,  increases  rapidly  with  A  when 
1  <  A  <  20  and  then  remains  high  for  20  <  A  <  100.  This 
means  that  the  modes  with  fast  increasing  ice  velocity  have 
short  band  spacing  (2  km  <  L/k  <  10  km). 

2.  For  fixed  H,  (1.5  m).  the  e-folding  acceleration  rate  o, 
increases  with  DTo,  since  the  larger  DTq  is.  the  stronger  the 
forcing. 


4' 
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Fig  4.  The  Ath  mode  of  ice  velocity. 
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Fig.  6-  The  e,lkyk  diagram  for  DTg  =  lO'C  and  (<i)  *  0.5  m,  (b) 

=  I  m,  (c)  H,  =  1.5  m.  (d)  //,-=■  2  m.  and  (e)  H,  =  2.5  m. 


3.  For  fi.xcd  DTf,  (lO'^C),  the  e-folding  acceleration  rate  a^ 
increases  with  a  decrease  in  Hf.  The  smaller  H,  is,  the  more 
unstable  the  coupled  system. 

4.  The  e-folding  acceleration  rate  for  the  most  unstable 
modes  is  around  1-4  x  10*’  s"‘.  and  the  doubling  time 
during  which  the  ice  doubles  its  speed  is 

T  =  In  2/<t,  (36) 

Therefore  the  doubling  time  for  most  unstable  modes  is 
around  5-20  hours. 

6.  Conclusions 

1.  This  ice-air  coupled  model  is  intended  to  depia  only 
the  mesoscalc  processes  of  ice-air  interaction  in  the  MIZ.  The 
synoptic  scale  pressure  gradient  may  additionally  produce  sur¬ 
face  winds  in  the  MIZ,  and  large-scale  ocean  currents  may 
drive  ice  drift.  These  processes  are.  however,  beyond  the  scope 
of  this  paper.  Nevertheless,  where  the  icc-to-opcn-watcr  tem¬ 
perature  gradient  is  strong,  the  mesoscalc  feedback  mechanism 
discussed  here  may  become  as  strong  as.  or  stronger  than,  the 
synoptic  scale  and  oceanic  forcings. 

2.  The  air-ice  feedback  mechanism  in  the  MIZ  generates 
the  alternation  pattern  of  ice  convergence  and  divergence 
zones.  Ice  floes  arc  thus  transported  from  the  divergence  area 


to  the  convergence  zone  and  are  eventually  rearranged  into  a 
band  struaurc  with  spacing  as  the  distance  of  two  adjacent 
convergence  zones. 

3.  The  alternation  pattern  of  ice  divergence/convergcnoc 
has  a  time  variation.  The  e-folding  acceleration  rate  of  ice 
velocity  depends  on  the  band  spacing  The  modes  with  a  fast 
increasing  ice  velocity  have  short  band  spacing  (2 
km  <  L/k  <  10  km),  which  agrees  well  with  the  observed  ice 
edge  band  spacing  [Hakkinen,  1986:  Martin  et  al..  1983: 
Muench  el  al..  1983;  Wadhams.  1983X 
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